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Abstract. We describe nonlinear quantum atom-light interfaces and nonlinear 
quantum metrology in the collective continuous variable formalism. We develop a 
nonlinear effective Hamiltonian in terms of spin and polarization collective variables 
and show that model Hamiltonians of interest for nonlinear quantum metrology can be 
produced in 87 Rb ensembles. With these Hamiltonians, metrologically relevant atomic 
properties, e.g. the collective spin, can be measured better than the "Heisenberg limit" 
oc 1/N. In contrast to other proposed nonlinear metrology systems, the atom-light 
interface allows both linear and non-linear estimation of the same atomic quantities. 
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1. Introduction 

In quantum metrology, a quantum state is prepared, evolves under the action of a 
Hamiltonian containing a parameter x of interest, and is measured. The parameter 
is estimated from the measurement outcome and knowledge of the system and 
Hamiltonian. In most problems, the Hamiltonian is assumed to act in the same way 
on each of N systems (e.g. atoms), and precision scales as Sx oc iV -1 / 2 for product 
states, and down to 5x oc iV -1 for entangled states ("Heisenberg limit" scaling) |JJ. 
A number of studies [2], [31 0] have considered also nonlinear quantum metrology, in 
which the Hamiltonian describes a fc-system coupling with strength x' . Remarkably, the 
scaling is Sx' oc N~ k+1 ^ 2 or Sx' oc N~ k for independent or entangled states, respectively 
[2]. Because this improves upon the best possible scaling for the linear case, it has been 
called "Super-Heisenberg" (SH) scaling [3]. Proposed implementations include scattering 
in Bose condensates p[| , Duffing nonlinearity in nano- mechanical resonators [5] , two-pass 
effective nonlinearity with an atomic ensemble |6], and Kerr nonlinearities [7J EJ |9]. 

Here we describe nonlinear metrology applied to measurement of collective 
spin variables of atomic ensembles. Atomic ensembles with long coherence-time 
internal degrees of freedom, e.g. nuclear spin, are essential elements of many 
quantum information and quantum metrology protocols including quantum memory 
fl()j. quantum non-demolition measurement [11] . spin squeezing [T2| fT3| Q3], and 
magnetometry |15j . 

We follow the approach of collective continuous variables (CCV), in which both 
light and atoms are described by macroscopic quantum variables. In the case of spinor 
atoms interacting with polarized light, the atoms are described by the collective spin 
F = f where f W is the spin of the z-th atom. F obeys the commutation relations 
[F x , F y ] = ihF z and cyclic permutations, and can itself be considered a macroscopic 
spin variable. The light is described by its electric field E = £ + £*, where £ is the 
positive-frequency part. The Stokes vector S with components Si = (£+, £l)<7i(£ + , £_) T 
where the subscript indicates plus/minus circular polarization, G{ are the Pauli matrices 
and o"o is the identity. As described by several authors [TBI [T71 [18], the electric dipole 
interaction h int = — E • d, taken in second order perturbation theory, gives rise to an 
effective (single-atom) Hamiltonian of the form 
^ £* -dJ0 fc ) (0 fc |d t -£ 
H<M= \ W k (1) 

= £*-a-£ (2) 

where d-f, dj. are the parts of the dipole operator causing upward/downward transitions, 
Sk is the detuning from resonance of the fcth state, and a is the tensor polarizability 
operator. The effective Hamiltonian for the ensemble H c g = ^ h^g can then be 
decomposed into irreducible tensor components as 

Hj^ = a 1 S z J z + a 2 {S X J X + S y J y ) , (3) 

plus terms in Sq which do not interact with the optical polarization. Here a.\ 2 describe 
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the vectorial, and tensorial components of the interaction, respectively, and the atomic 
collective variable is J = X1J W with 3x = (//-//) A j v = (f x f v + fyfx) A 
j 2 = / z /2 and j = A The ratio of the ati can be tuned by adjusting the optical 
frequency u, giving a variety of Hamiltonians interesting for quantum information tasks 

To apply this formalism to nonlinear metrology, we generalize the CCV method 
to the nonlinear optics regime, i.e., we include higher-order processes in the effective 
Hamiltonian. For this purpose, naive application of higher-order perturbation theory 
fails due to the appearance of vanishing resonance denominators, and degenerate 
perturbation theory [19] is required. We present the method by way of an example, 
the D 2 line of 87 Rb, one of the most used transitions for atom-light interactions. 

2. Derivation of the Effective Hamiltonian 

We consider the 5 2 Si/ 2 — > 5 2 P 3 / 2 transition (the D 2 transition at 780 nm). The ground 
states are \F,m F ) with (F,m F ) = (1, -1), (1, 0), (1, 1), (2, -2), . . . , (2, 2). The excited 
states are \F', irtpi) with (F', mpi) = (0, 0), (1, —1), (1,0),..., (3, 3). We use these states 
as a basis, with the ground states preceding the excited states. We calculate the single- 
atom Hamiltonian h c g, and note that the ensemble Hamiltonian H e ^ = Y2i ^is * s f° un d 
simply by replacing single-atom operators such as j with collective operators J = jW. 
The unperturbed Hamiltonian is ho = h^2 l U)i \l) or in matrix notation 

h = h(u F=1 I 3 © 0J F=2 h © oj f > =q Ii © oj f > =1 J 3 

®U) F , =2 h ®uj f , =3 I 7 ). (4) 

where © indicates a direct sum, and is the identity matrix of dimension d. We choose 
the origin of energy such that u F= i = 0, and define A = bJ F=2 - We work in a frame 
rotating with the laser frequency u = u) F >=o + 5. In this frame, the Hamiltonian is 

ho = h(0I 3 © A/ 5 © Soh © <Ji/ 3 © S 2 h © hh) (5) 

where 5 pi = u F i — oo. 

In the rotating wave approximation, the single-atom perturbation v = h int = — E d 
is approximated as v ~ £ d^+£* d^. If S± are the amplitudes for the sigma-plus/minus 
components, respectively, of £ , then 

(F',m F i\V\F,m F ) = £ q (F' , m F i \ er q \F, m F ) (6) 

with q = m F i — mp. Note that q = transitions (7r-transitions) are not considered 
because the 2-propagating beam cannot contain this polarization. The dipole matrix 
elements are related to the "matrix element" (J| \er q \ \ J') = Djj> w 3.58410~ 29 C • m by 
angular- momentum addition rules. We follow the conventions given in Steck |20j. In 
this way, we arrive to the perturbation Hamiltonian 
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where 0^ = 0/^ and = \^5Djji x 
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To obtain the effective Hamiltonian, we follow Klein |19| . The notation of that 
work is somewhat obscure, so for ease of understanding we repeat the main results. 
From Equation (A7) of that work, we have the t-order contribution to the effective 
Hamiltonian 



i {t) 



Where k±, 



( 9 ) 
{k} 

. , k t -i are non- negative integers, the A are real coefficients, the O, denoted 
kt)" by Klein, are operators, and the sum is taken over all {k} satisfying 
— 1. The A are given in Table I of that work and the O are given in Equation 



(Al) as 

O ikl ,..., h _ l} = P VR {k ^VR {k2) . . . VR^VPo (10) 
with P being the projector onto the degenerate subspace and by Equation (II. A. 5) 

P k = 



E -H(°) _ 

where E Q is the energy of the degenerate subspace. In our case we have chosen E = 0. 

We can then directly calculate the second- and fourth-order contributions. We are 
only concerned with h e g as it acts on the F = 1 subspace, that is, with a 3 x 3 matrix, 
and it is convenient to express it in terms of the pseudo-spin components jo, jx, jy, jz 
and the Stokes components So, S x , S y , S z defined above. Summing the second-order 




Figure 1. (color online). Spectra of the terms of the effective Hamiltonian. First 
two curves from the top, left axis: continuous, o^ 1 ); dashed, a^ 2 \ Lower curves, right 
axis: continuous. /3 (1) ; dotted, 0$; dashed, /3 (2) ; dot-dashed, /3 ( 7 0) . Detuning [MHz] is 
relative to the transition F = 1 —> F' = of 87 Rb D2 transition. Points A, B, indicate 
detunings at which or /?W vanish. 



contributions we find of equation [3] With B = — Djj,/485o5i52h, 

a® = 5(55 5i - 55 5 2 - 4^5 2 ) (12) 
a (2) = 5(5 5i - 55 5 2 + 45i5 2 ) (13) 
Similarly, the fourth-order contribution is, dropping terms in Sq, 

H% = Pf&zh + PfslN A + pVSoSzJz 

+ P ( - 2) S (S X J X + S Y J Y ). (14) 

Note that the term in N A arises because contains a self-rotation term of the form 
/3^1 S'|P fn= o where P m =o is a projector onto the state \F — 1, tuf — 0). We express this 
in terms of J and N A using £\ P%L = J2i(4 ~ 3o ) = N A - J . 

With C = Djj//11525q<5^5|A?i 3 , the coefficients, shown graphically in Figure [U are 

pf = C (\25l8l5l - ±815x81 + 12<Jg<J?A - 10<5 3 5 2 <5 2 A - 12<5 2 5 3 <5 2 A 

-lOc^M^A - l25 G 5l5 2 2 A + 20c^M 2 A + 205 5 2 5 2 A) (15) 

pf = C (-125 3 <5^ 2 - 2A5l5l5l + 45 3 M 2 3 ) (16) 
= C (-95l5l5 2 + 6<5 3 5 2 <5 2 + 3<5g5i<5| + 355 3 <5 3 A - 5<5 3 <5 2 5 2 A 
-4<5 2 5 3 5 2 A - 5<Jj}£i<Sf A - 45 5 3 5 2 A - 25^ A - 20<^itff A 
-205 5 2 ^A - 16(5?<J|A) (17) 
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+16<5g<y?<5 , 2 A - 15<$5i5|A + 165 5^2 A 
+16<??<JfA) 



25<^A 



3. Application to nonlinear metrology 

The (3 terms are nonlinear in S, indicating a photon-photon interaction. We expect 
these terms to describe polarization effects of fast electronic nonlinearities including 
saturation and four-wave mixing. As in the linear case, the frequency dependence of the 
/3 terms provides considerable flexibility in designing a light-matter interaction. Applied 
to quantum metrology, these terms produce SH scaling, because they are nonlinear in the 
S collective variables, while the atomic variables Jj, play the role of the parameter. 
The /?(°) and (3^ terms are analogous to Hamiltonians considered by Boixo et al. [2]. 
The (3^ term oc SqS z , in particular, achieves SH scaling without input or generated 
entanglement The (3^ terms describes a nonlinear tensorial contribution, and does 
not appear to have been considered yet for nonlinear metrology. 



4. Quantum Noise 

To understand the quantum noise in this system, we define polarization operators 
Si = \{a} + , a^)<7j(a + , a_) T , where a± are annihilation operators for the ± circular 
polarizations of a mode defined by the pulse shape. These obey angular momentum 
commutation relations [Si,Sj] = iEijkSk and are related to the Stokes parameters by 
Si = Si/27, where 7 = HuZq/2TA is the single-photon intensity, T is the pulse duration, 
A is the beam area, and Z is the impedance of free space. The total number of photons 
is N L = 2S . For a typical input, a coherent state, ( (Sx, Sy, Sz) ) = (So, 0,0) and 
var(Si) = S /2. 

Evolution under this effective Hamiltonian produces, to first order in the interaction 
time r, 

g(o«t) = £H + ^ a (i) + pW^jSP J| n) (19) 

plus terms containing S^J x n ^ that are negligible for the given input coherent state of 
the light. This evolution physically corresponds to a paramagnetic Faraday rotation 
of the input linear polarization. In a metrological scheme one would measure this 
polarization rotation and from it estimate the atomic variable Jz- 

For small rotation, i.e. = S^ ut) /S x n) = 7T/T 1 (a (1) + /3 (1) 7S )j| n) < 1, we note 
that the input polarization noise dominates: var(Sy Ut ^) = var(Sy 11 ' ) ) + </> 2 var(Sjf ) ~ 
var(Sy n ^) = So/2, and that the signal-to-noise ratio equals one when ( Sy Ut ^ ) 2 = 
var(Sy^), i.e., when 

^S 2 (« (1) +/3 (1) 7S ) 2 (j| n) ) 2 = |. (20) 
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We can identify the value of that solves Eq. ( 120]) as the sensitivity, or precision 
of the estimation, SJz- We find 

5J Z = hlrjia^N 1 ^ + ^l N 3 L /2 )\-\ (21) 

Thus the sensitivity will have a transition from shot-noise to SH scaling with increasing 
Nl- As indicated in Figure [1], there are points in the spectrum where either or (3^ 
vanish, allowing pure nonlinear or pure linear estimation of the same atomic variable. 

In another scenario, an unpolarized input state ( (£>x, Sy, Sz) ) = (0, 0, 0) gives rise 
to dynamics dominated by the (3^ terms oc S\ , sometimes called the "one-axis twisting 
Hamiltonian." This describes a self-rotation of the optical polarization, and can be used 
to generate polarization squeezing and also to obtain sensitivity scaling as N^ 3 ^ 2 in the 
estimation of pf^Jo + (3^Na, using an entanglement-generating strategy described in 
reference [3|. 

5. Conclusion 

We have generalized the formalism of continuous collective variables to the nonlinear 
regime. The resulting nonlinear effective Hamiltonian includes several distinct nonlinear 
couplings with strengths widely tunable via the probe light frequency. This allows the 
production of model Hamiltonians proposed for nonlinear metrology, including both 
models that generate entanglement and those which achieve super-Heisenberg scaling 
without entanglement. Similar nonlinear probing techniques could improve optical 
probing of atomic clocks |21 | IT3 ] and atomic magnetometers [221 IT5| 123 ] . Unlike previous 
proposals, the atomic ensemble system allows both linear and nonlinear estimation of 
the same atomic variables. 
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